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We report on the implementation of a time-multiplexed click detection scheme to probe quantum 
correlations between different spatial optical modes. We demonstrate that such measurement setups 
can uncover nonclassical correlations in multimode light fields even if the single mode reductions are 
purely classical. The nonclassical character of correlated photon pairs, generated by a parametric 
down-conversion, is immediately measurable employing the theory of click counting instead of low- 
intensity approximations with photoelectric detection models. The analysis is based on second- and 
higher-order moments, which are directly retrieved from the measured click statistics, for relatively 
high mean photon numbers. No data postprocessing is required to demonstrate the effects of interest 
with high significance, despite low efficiencies and experimental imperfections. Our approach shows 
that such novel detection schemes are a reliable and robust way to characterize quantum-correlated 
light fields for practical applications in quantum communications. 

PACS numbers: 42.50.-p, 42.50.Ar, 03.65.Wj 


Introduction.- Certifying quantum features of light 
is one key requirement for optical implementations of 
quantum-information technology mm- This requires, 
on the one hand, reliable sources for correlated quan¬ 
tum light and, on the other hand, appropriate detection 
schemes mm- Since correlations between different de¬ 
grees of freedom can have different origins in quantum 
optics, they might be covered by classical statistical op¬ 
tics, or they are genuine quantum properties having no 
such classical counterpart. Using the classical theory of 
coherence, one way to discern quantum from classical ef¬ 
fects has been introduced independently by Glauber 
and Sudarshan [5]. 

A quantum-state characterization is often based on the 
photon number distribution, see, e.g.. Refs. Clll]. The 
corresponding photoelectric detection theory yields Pois- 
sonian statistics for coherent light. However, detectors 
that directly measure photon numbers are typically not 
available or require advanced data postprocessing liiin]. 
Today, quantum states with low mean photon number 
are often detected with avalanche photodiodes (APDs) 
in the Geiger mode, which basically produce a “click” if 
any number of photons is absorbed, and remain silent 
otherwise. A uniform splitting of a radiation field with 
many photons into portions of lower intensities, each be¬ 
ing measured with an APD, can extent the knowledge 
about the signal, for example, to discriminate single- 
and two-photon events. Various kinds of such photon- 
number-resolving detectors have been implemented to 
demonstrate nonclassical features of radiation fields, e.g., 
in Refs. EHn], or for characterizing the non-Poissonian 
behavior of the click statistics [181 [E] ■ One particular 
realization of such a scheme that requires only a small 
number of optical elements is so-called time-bin multi- 
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FIG. 1: (Color online) The realization of two TMDs in a 
single device is depicted. One part (upper, green input) of 
the correlated signal field is delayed in an optical fiber before 
entering the device. Within the TMD, any signal is split by 
a 50/50 beam splitter followed by another delay line in one 
output. This multiplexing is done three times yielding 2^ 
separated time bins for each (upper green and lower blue) 
input field which are measured with two APDs. 


plexing detectors (TMDs) [Ill[20l|2T], cf. Fig. Based 
on these TMDs, one can reconstruct nonclassical features 
of quantum light fields or higher-order correlation 

functions [24] . 

The proper theoretical detection model for such de¬ 
vices is a quantum version of the binomial statistics [25j . 
It was also shown that approximating these statistics 
with a Poissonian distribution and applying quantum¬ 
ness probes of the photon statistics can yield fake signa¬ 
tures of nonclassicality for classical fields - even if the 
number of registered signal photons is one order of mag¬ 
nitude below the number of time bins. To eliminate these 
errors, nonclassicality criteria in terms of moments of 
the click-counting statistics have been proposed to di¬ 
rectly uncover quantum light with click counters without 
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data postprocessing [^. For example, the notion of sub¬ 
binomial light has been introduced |27) and experimen¬ 
tally demonstrated [28l [29] . 

In the present work, we report on the characteriza¬ 
tion of a bipartite quantum-correlated light source using 
click detectors only. We demonstrate that click detec¬ 
tion is capable of directly verifying nonclassicality even 
for imperfect experimental settings. For our parametric 
down-conversion (PDC) based light source, this approach 
correctly shows classical single-mode correlations and, at 
the same time, it uncovers quantum correlations between 
the modes. We show that our simple treatment to infer 
these quantum features works for a broad range of pump 
powers. In particular, it works increasingly well for in¬ 
creasing pump powers, where the often used weak signal 
approximations with the photoelectric detection theory 
completely break down, as we will show later. 

Nonclassical moments of the click statistics.- The 
probability to measure /ca clicks within the N = 8 time 
bins assigned to the signal A together with clicks from 
the signal B, cf. Fig. is described through the joint 
click-counting statistics [IHl HB] : 

CkA.kB 

^ ^ ( 1 ) 

with : • : denoting the normally ordering prescription and 
TTii = for the modes i = A, B. In general, F 

can be an unknown detector response, being a function 
of the photon number operators hi. For example, a linear 
form of the response function is r{hi/N) = rjhi/N + i/, 
with rj and v being the quantum efficiency and the dark 
count rate, respectively. 

In Ref. |2B| it has been demonstrated that the matrix of 
click moments is non-negative for any classical 

light field, 

0 < ( 2 ) 


with SiAi = ■ ,Kil2 < N/2 for even Ki and N. The 

superscript {Ka, Kb) defines the highest moment of each 
subsystem within the matrix M, see also Supplement. 
For instance, the single-mode and bipartite, second-order 
matrices of click moments are 


M( 2 . 0 ) = 


and = 
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respectively. The needed moments can be directly re¬ 
trieved from the measured click-counting statistics [ 2 Bj, 
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One way to probe the character of nonclassical correla¬ 
tions, i.e., violating inequality ([^, can be done as follows. 
We have nonclassical ifth-order click correlations if 

> 0, > 0 , and ^ q. (5) 

This means that both RTth-order single-mode marginals 
are classical and the bimodal, RTth-order correlation ma¬ 
trix is nonclassical, i.e., it has at least one negative eigen¬ 
value. In order to genuinely certify such nonclassical cor¬ 
relations, it is sufficient to consider the minimal eigenval¬ 
ues of the click-moment matrices. Say /a, /b, and /as 
are the normalized eigenvectors to the minimal eigenval¬ 
ues CA, es, and cab of ^nd ^ 

respectively. Now, definition ([^ is rewritten as 

eA = flM<^^’°^fA > 0, CB = > 0, 

and CAB = fAB^^^'^^fAB < 0. 

This method will serve as our approach to determine 
i^rth-order quantum correlations between the subsystems 
A and B, see Supplement for further details. 

Implementation and model.- The states under study 
are produced in a type II PDC in a periodically poled, 
8 mm long KTP waveguide. The PDC process is pumped 
with 1 ps long pulses at a repetition rate of 70 kHz and 
a wavelength of 768 nm coming from a Ti:Sapphire laser. 
The states are generated in two orthogonally polarized 
signal and idler modes at 1536 nm. In the regime of less 
than one photon per pulse, the source has been charac¬ 
terized in Ref. |30j with emphasis on the modal proper¬ 
ties, and genuine single-mode operation of the PDC has 
been shown. Behind the waveguide, broadband spectral 
filters are used to suppress the pump and unwanted back¬ 
ground outside of the PDC spectral region. Signal and 
idler modes are then split at a polarizing beam splitter 
and coupled into a two-mode TMD consisting of a fiber 
network with 50/50 beam splitters and a pair of InGaAs 
APDs as depicted in Fig. For each pump setting, we 
record the full time tagging data for 10 min and trans¬ 
late them into the click statistics of all possible 8 x 8 
click events. Experimental imperfections are unbalanced 
beam splitters within our TMD, i.e., splitting ratios that 
slightly deviate from 50/50, or after-pulsing effects [5T] . 
We minimized both effects through a careful setup prepa¬ 
ration. 

The PDC may be formulated in terms of the effective 
Hamiltonian H = iK'jafw + H.c., where k is a coupling 
constant, 7 is the coherent amplitude of the pump beam, 
and a^( 6 'f) is the creation operator of the mode A[B). 
An idealized, perfect unitary evolution for this process 
yields the two-mode squeezed-vacuum state, 

Ij) = (7) 

where ^ > 0 is proportional to the square root of the 
pump power P, since oc 7 and P oc 7 ^. The noise 
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FIG. 2: (Color online) The joint click-counting statis¬ 
tics 0 for the state 0. ^=1, and a perfect linear response, 
r{h/N)=h/N . The probability for no clicks is cut, co,o=0.42. 
Even though we have perfect photon correlations, the click¬ 
counting distribution includes off-diagonal terms. 


suppression of the squeezed quadrature X is ([AX]^) = 
e“^^([AX]^)vac [Z] and Supplement. Because of the pair¬ 
wise generation of photons, the state 0 has perfect 
photon-number correlations. This means that whenever 
a certain number of photons is present in one mode, the 
same number of photons occurs in the other mode. How¬ 
ever, the click-counting statistics include off-diagonal el¬ 
ements CkA,kB 7^ 0 foi' 7 ^ ks, even for a perfect detec¬ 
tion without dark counts and unit efficiency, cf. Fig. 
This difference between click-counting and photoelectric 
detection is due to a hnite probability that more than one 
photon can end in the same time bin. It is worth men¬ 
tioning that the single-mode reduced states trA|^)(^| and 
t''_B|C)(CI are classical thermal states, and that the total 
number of photons is h = (■Clh-Ti + nslO = 2sinh^ 

Second-order correlations.- First, let us focus on 
second-order click correlations, cf. Eq. 0, which include 
the information about the mean values, the variances, 
and the covariance of the joint click-counting statis¬ 
tics [5S]. In Fig. 0 we plot the measurement results. 
Using the approach in Eq. 0, the minimal eigenval¬ 
ues CAB (top), CA (bottom, left), and cb (bottom, right) 
are shown in their dependence on the energy per pulse, 
.Upump = H/ZOkHz. The energy can be experimentally 
controlled. The single-mode matrices are non-negative, 
CA > 0 and es > 0, whereas the cross correlations are 
nonclassical, cab < 0. Thus, we have verified the quan¬ 
tum nature of the second-order click correlations between 
the spatial modes A and B. 

In order to compare our measured results, a simple the¬ 
oretical model is used. We assume that the pure state 0 
is generated, and the detectors are described via a plain 
linear response function: T(h/N) = njh/N-{-u. As we dis¬ 
cussed above, the parameter characterizing the state 
1^), depends on the pump power P: ^ = ^oVP- The pro¬ 
portionality constant the quantum efficiency 77 , and 
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FIG. 3: The upper plot shows the minimal eigenvalue cab 
(black bars) of the experimentally obtained matrix of mo¬ 
ments depending on the pump energy Epu^p. Error 

bars are given as gray areas. The dashed curve is the theoret¬ 
ical prediction. The inset depicts the continuation of the the¬ 
oretical curve for higher energies including saturation effects. 
The lower plots show the minimal eigenvalues ca and es of 
(left) and (right), respectively, a ten standard 

deviations error bar, and the theoretical prediction. Since the 
bimodal correlations are negative and the single-mode reduc¬ 
tions are non-negative, we successfully determined nonclassi¬ 
cal correlations between the modes A and B. 


the dark count rate have been fitted using the standard 
method of least squares. They are r] = 9.6%, o = 0.51, 
and ^0 = 0.087 (/iW)“^/^. With these values, we es¬ 
timate the total mean photon numbers to be 0.9... 15 
photons for the pump powers 50... 403 or energies 
0.7... 5.8 nJ. Already this simplified model yields a good 
agreement with the measured results, which highlights 
the excellent performance of the engineered PDC source, 
cf. the dashed lines in Fig. 0 

The inset in the upper plot shows the extrapolation 
of the quantum correlations, cab < 0, for higher pump 
energies. At some point, the mean photon number is 
so high, that these correlations saturate and eventually 
vanish. A high squeezing level and a classical laser light 
with a large coherent amplitude result in the same sig¬ 
nal - all time bins are occupied with a large number of 
photons. Therefore, at high photon numbers, the signals 
of nonclassical and classical states cannot be discrimi¬ 
nated. Thus, the recorded quantum correlations of the 
former state must decrease due to the saturation of click¬ 
counting devices, which is automatically included in the 
click-counting theory. Imperfections are also studied in 
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the Supplement. 

Let us emphasize again that the states have been gen¬ 
erated for pump powers ranging over almost one order 
of magnitude. Verifying nonclassical photon-photon cor¬ 
relations in this comparably large domain is typically 
considered a challenging task, but can easily be accom¬ 
plished with our TMD click counters. In addition, since 
the method of nonclassical click moments is independent 
of the state, the verified quantum correlations can be cer¬ 
tified even if the pump power was completely unknown. 

Higher-order correlations.- Let us study higher-order 
quantum correlations. They become particularly mean¬ 
ingful when second-order criteria fail to properly charac¬ 
terize the state [32]. For example, the third- and fourth- 
order moment relate to the so-called skewness and the 
flatness (or kurtosis), respectively. The highest possible 
order of moments one can infer from N = 8 time bins 
per mode is IF = 8 in Eqs. ([^ and Q, which yields a 
full characterization of the click-counting statistics. 

The bound for a classical signal, > 0, is given 

by the eigenvalue Cd = 0 [55]. Thus, the signed distance, 
in units of standard deviations, of the experimentally ob¬ 
tained minimal eigenvalue e = e ± Ae to this classical 
bound leads to a signed significance. 


representing a signed relative error. A negative signifi¬ 
cance E < 0 verifies the nonclassicality with a significance 
of |E|. Typically, E < —3 is a significant verification of 
the negativity, whereas E « 0 cannot be distinguished 
from the classical bound 0 . 

In Fig. [^ the significance levels of the full eighth or¬ 
der quantum correlation within are given. The 

single-mode, signed significances for and 

can be found in the Supplement. There are no significant 
eighth order, single-mode correlations > 0 

- the largest single-mode negativities are of the order of 
E ^ —10“®. Additionally, the significance of negativi¬ 
ties in M(®’®\ cf. Fig. S' are in the range of 2.7... 10.6 
for the energies 0.7...5.8nJ, respectively. Hence, for 
most of the energies, a significant eighth-order nonclassi¬ 
cal correlation between the modes A and B is certified. 
Remarkably, the significance even increases with the en¬ 
ergy, which is due to an improved signal-to-noise ratio 
with increasing mean photon numbers, because the no¬ 
click event has a much lower probability in this regime. 

On the one hand, one would typically not use such 
comparably high intensities in our measurement setup, 
when analyzing the data with the inappropriate photo¬ 
electric detection model. In this case, the single-mode 
signed significances evaluate to — 2.4... —13, cf. Supple¬ 
ment, suggesting fake nonclassicality which worsens with 
increasing pump power. On the other hand, our consis¬ 
tent treatment in terms of the click statistics correctly 
identifies higher-order bimodal correlations while show- 



FIG. 4: The signed significances E of the largest negativities 
of the full matrix of click moments are shown. The 

highest significance levels of nonclassicality between the sub¬ 
systems are verified for our highest energy. In this region the 
impact of the click-counting theory is most pronounced. 


ing, as expected for our source, no nonclassicality in the 
single-mode marginals. 

Here, we studied photon-photon correlations in terms 
of click-counting moments for high mean photon num¬ 
bers. Even though we have an incomplete knowledge 
about the state - no phase information and only a fi¬ 
nite number of bins - we are able to uncover quan¬ 
tum correlations. The lack of phase resolution compares 
to a fully phase randomized version of the state 0 . 
which has been shown to include nonclassical correla¬ 
tions without entanglement |33j . being certified with the 
present method. The measurement scheme may be com¬ 
plemented for phase-sensitive measurements |34j , to infer 
squeezing or entanglement. 

Conclusions. - We set up a click correlation measure¬ 
ment scheme using time-multiplexed detectors for prob¬ 
ing nonclassical correlations of a bipartite quantum light 
field. We could directly infer nonclassical correlations 
from the measured joint click-counting statistics with 
high significances. This was possible despite a low de¬ 
tection efficiency, estimated as 10 %, and over one order 
of magnitude of intensities. Note that the bare click data 
have been analyzed without corrections for imperfections 
or related postprocessing techniques. We compared our 
results with a theoretical model and obtained a very good 
agreement. In particular, the expectation of this model 
- all orders of single-mode correlations do not exhibit 
nonclassicality, but the two-mode correlations do - is 
correctly demonstrated with our method. This under¬ 
lines the functionality of our device for verifying genuine 
quantum correlations between spatial optical modes. 

Our proposed measurement together with the click- 
detection analysis is a robust and efficient tool to char¬ 
acterize quantum light. We believe that this simplicity 
of click detection - being solely a collection of probabili¬ 
ties of coincidence clicks - paves the way towards real-live 
implementations of quantum communication protocols in 
optical systems. The present approach may be further 
generalized to handle more complex types of quantum- 
correlated, multimode radiation fields. 








5 


Acknowledgements.- J.S., M.B., and W.V. acknowl¬ 
edge financial support by Deutsche Forschungsgemein- 
schaft through SFB 652 and VO 501/22-1. G.H., B.B., 
V.A., and C.S. acknowledge funding from the European 
Community’s Seventh Framework Program FP7/2001- 
2013 under Grant Agreement No. 248095 through the 
Integrated Project Q-ESSENCE. 

Note added.- We recently became aware of a related 
paper in preparation by the group of I. A. Walmsley and 
co-workers [55] . 


[1] P. Kok, W. J. Munro, K. Nemoto, T. C. Ralph, J. P. 
Dowling, and G. J. Milburn, “Linear optical quantnm 
compnting with photonic qubits”. Rev. Mod. Phys. T9, 
135 (2007). 

[2] A. Datta, L. Zhang, N. Thomas-Peter, U. Dorner, B. 
J. Smith, and I. A. Walmsley, “Quantum metrology with 
imperfect states and detectors”, Phys. Rev. A 83, 063836 
( 2011 ). 

[3] G. S. Buller and R. J. Collins, “Single-photon generation 
and detection”, Meas. Sci. Technol. 21, 012002 (2010). 

[4] M. D. Eisaman, J. Fan, A. Migdall, and S. V. Polyakov, 
“Single-photon sources and detectors”. Rev. Sci. In¬ 
strum. 82, 071101 (2011). 

[5] R. J. Glauber, “Photon Correlations”, Phys. Rev. Lett. 
10, 84 (1963); R. J. Glauber, “Coherent and Incoherent 
States of the Radiation Field”, Phys. Rev. 131, 2766 
(1963). 

[6] E. C. G. Sudarshan, “Equivalence of Semiclassical and 
Quantum Mechanical Descriptions of Statistical Light 
Beams”, Phys. Rev. Lett. 10, 277 (1963). 

[7] G. S. Agarwal, Quantum Optics (Cambridge University 
Press, Cambridge, 2013). 

[8] L. Mandel and E. Wolf, Optical Coherence and Quantum 
Opties (Cambridge University Press, Cambridge, 1995). 

[9] O. Haderka, J. Pefina, Jr., M. Hamar, and J. Pefina, “Di¬ 
rect measurement and reconstruction of nonclassical fea¬ 
tures of twin beams generated in spontaneous parametric 
down-conversion”, Phys. Rev. A 71, 033815 (2005). 

[10] V. N. Starkov, A. A. Semenov, and H. V. Gomonay, “Nu¬ 
merical reconstruction of photon-number statistics from 
photocounting statistics: Regularization of an ill-posed 
problem”, Phys. Rev. A 80, 013813 (2009). 

[11] E. Waks, E. Diamanti, B. C. Sanders, S. D. Bartlett, 
and Y. Yamamoto, “Direct Observation of Nonclassi¬ 
cal Photon Statistics in Parametric Down-Conversion”, 
Phys. Rev. Lett. 92, 113602 (2004). 

[12] A. Allevi, M. Bondani, and A. Andreoni, “Photon- 
number correlations by photon-number resolving detec¬ 
tors”, Opt. Lett. 35, 1707 (2010). 

[13] G. Puentes, A. Datta, A. Feito, J. Eisert, M. B. Plenio 
and 1. A. Walmsley, “Entanglement quantification from 
incomplete measurements: applications using photon- 
number-resolving weak homodyne detectors”. New J. 
Phys. 12, 033042 (2010). 

[14] J. F. Dynes, Z. L. Yuan, A. W. Sharpe, O. Thomas, 
and A. J. Shields, “Probing higher order correlations of 
the photon field with photon number resolving avalanche 
photodiodes”. Opt. Express 19, 13268 (2011). 


[15] P.-A. Moreau, J. Mougin-Sisini, F. Devaux, and E. Lantz, 
“Realization of the purely spatial Einstein-Podolsky- 
Rosen paradox in full-field images of spontaneous para¬ 
metric down-conversion”, Phys. Rev. A 86, OIOIOI(R) 
( 2012 ). 

[16] L. Dovrat, M. Bakstein, D. Istrati, E. Megidish, A. 
Halevy, L. Cohen, and H. S. Eisenberg, “Direct observa¬ 
tion of the degree of correlations using photon-number- 
resolving detectors”, Phys. Rev. A 87, 053813 (2013). 

[17] G. Donati, T. J. Bartley, X.-M. Jin, M.-D. Vidrighin, 
A. Datta, M. Barbieri, and 1. A. Walmsley, “Observing 
optical coherence across Pock layers with weak-field ho¬ 
modyne detectors”, Nat. Commun. 5, 5584 (2014). 

[18] M. J. Fitch, B. C. Jacobs, T. B. Pittman, and J. D. Fran- 
son, “Photon-number resolution using time-multiplexed 
single-photon detectors”, Phys. Rev. A 68, 043814 
(2003). 

[19] J.-L. Blanchet, F. Devaux, L. Furfaro, and E. Lantz, 
“Measurement of Sub-Shot-Noise Correlations of Spa¬ 
tial Fluctuations in the Photon-Counting Regime”, Phys. 
Rev. Lett. 101, 233604 (2008). 

[20] J. Rehacek, Z. Hradil, O. Haderka, J. Pefina Jr., and M. 
Hamar, “Multiple-photon resolving fiber-loop detector”, 
Phys. Rev. A 67, 061801 (R) (2003). 

[21] D. Achilles, C. Silberhorn, C. Sliwa, K. Banaszek, and 
1. A. Walmsley, “Fiber-assisted detection with photon 
number resolution”. Opt. Lett. 28, 2387 (2003). 

[22] K. Laiho, K. N. Cassemiro, D. Gross, and C. Silberhorn, 
“Probing the Negative Wigner Function of a Pulsed Sin¬ 
gle Photon Point by Point”, Phys. Rev. Lett. 105, 253603 
( 2010 ). 

[23] G. Harder, C. Silberhorn, J. Rehacek, Z. Hradil, L. 
Motka, B. Stoklasa, and L. L. Sanchez-Soto “Time- 
multiplexed measurements of nonclassical light at tele¬ 
com wavelengths”, Phys. Rev. A 90, 042105 (2014). 

[24] M. Avenhaus, K. Laiho, M. V. Chekhova, and C. Silber¬ 
horn, “Accessing Higher Order Correlations in Quantum 
Optical States by Time Multiplexing”, Phys. Rev. Lett. 
104, 063602 (2010). 

[25] J. Sperling, W. Vogel, and G. S. Agarwal, “True pho¬ 
tocounting statistics of multiple on-off detectors”, Phys. 
Rev. A 85, 023820 (2012). 

[26] J. Sperling, W. Vogel, and G. S. Agarwal, “Correlation 
measurements with on-off detectors”, Phys. Rev. A 88, 
043821 (2013). 

[27] J. Sperling, W. Vogel, and G. S. Agarwal, “Sub-Binomial 
Light”, Phys. Rev. Lett. 109, 093601 (2012). 

[28] T. J. Bartley, G. Donati, X.-M. Jin, A. Datta, M. Bar¬ 
bieri, and I. A. Walmsley, “Direct Observation of Sub- 
Binomial Light”, Phys. Rev. Lett. 110, 173602 (2013). 

[29] R. Heilmann, J. Sperling, A. Perez-Leija, M. Grafe, M. 
Heinrich, S. Nolle, W. Vogel, and A. Szameit, “Divide- 
and-Conquer: An integrated photon-counting scheme”, 
arXiv: 1502.04932 [quant-ph]. 

[30] G. Harder, V. Ansari, B. Brecht, T. Dirmeier, C. Mar- 
quardt, and C. Silberhorn, “An optimized photon pair 
source for quantum circuits”. Opt. Express 21, 13975 
(2013). 

[31] Y. Kang, H. X. Lu, and Y.-H. Lo “Dark count proba¬ 
bility and quantum efficiency of avalanche photodiodes 
for single-photon detection”, Appl. Phys. Lett. 83, 2955 
(2003). 

[32] G. S. Agarwal and K. Tara, “Nonclassical character of 



6 


states exhibiting no squeezing or sub-Poissonian statis¬ 
tics”, Phys. Rev. A 46, 485 (1992). 

[33] A. Ferraro and M. G. A. Paris, “Nonclassicality Crite¬ 
ria from Phase-Space Representations and Information- 
Theoretical Constraints Are Maximally Inequivalent”, 
Phys. Rev. Lett. 108, 260403 (2012); E. Agndelo, J. Sper¬ 
ling, and W. Vogel, “Quasiprobabilities for multipartite 
quantum correlations of light”, Phys. Rev. A 87, 033811 
(2013). 


[34] J. Sperling, W. Vogel, and G. S. Agarwal, “Balanced ho¬ 
modyne detection with on-off detector systems: Observ¬ 
able nonclassicality criteria”, Europhys. Lett. 109, 34001 
(2015); A. Luis, J. Sperling, and W. Vogel “Nonclassi¬ 
cality Phase-Space Functions: More Insight with Fewer 
Detectors”, Phys. Rev. Lett. 114, 103602 (2015). 

[35] T. Bartley et al. (unpublished) 



Supplemental Material 

Uncovering Quantum Correlations with Time Multiplexed Click Detection 


J. Sperling, M. Bohmann, and W. Vogel 

Arbeitsgruppe Theoretische Quantenoptik, Institut fiir Physik, Universitat Rostock, D-18051 Rostock, Germany 


G. Harder, B. Brecht, V. Ansari, and C. Silberhorn 
Integrated Quantum Optics Group, Applied Physics, 

University of Paderborn, 33098 Paderbom, Germany 

The supplement is structured as follows. In Secs. I and II, we discuss how to infer the moments 
and error bars of the measured joint click-counting statistics. In Sec. Ill, we give the details on 
the determination of classical and nonclassical matrices of moments in terms of eigenvalues and 
their significances. We present all the results of our analysis in Sec. IV. Fake quantum effects due to 
photoelectric counting approximations are discussed in Sec. V. In Sec. VI, we analyze the theoretical 
model together with a discussion of imperfections. 


I. MEASURING THE CLICK STATISTICS 


All we need is a plain device producing a binary information referred to as “click” and “no-click” or “on/off” events. 
This can be formally written as the expectations value of a so-called detector response function r(h), were h is the 
photon number operator. Thus the positive operator-valued measure is given by the elements 

Hciick = i — and Hno-ciick = . (1) 


Splitting a signal into N = 8 modes of equal intensities - each of them measured with this binary “on/off” device - 
yields the click counting statistics of two detectors [1]: 


CkA,kB 




N-kB 


1 — e 


-TifiB/N) 


\ kB 

) --k 


( 2 ) 


where ki (i = A, B) denote the numbers of coincidence clicks from both detectors. Note that we assume that both 
detectors A and B have the same characteristics, as it is true for our implementation. A typical example is a linear 
response function, r(h/V) = rjfi/N + v, with r] being the quantum efficiency, and v the dark count rate. 

It is worth to stress that CkA,kB is the probability that we have kA clicks from the first detector system simultaneous 
with ks clicks of the second one. The experimentally recorded number of coincidence events (fc^, ks) may be labeled 
as CkA,kB- Hence, we get the relative frequencies 
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^kA,kB 
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with C = E E 

k^—O kB—0 


(3) 


denoting the total number of events. Later on, the single mode marginals will have some importance. They are given 


by c' 


kA 


_ ^-^iV 
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0 ^kA ikE 


1 exp 

and c. P = Y.kA=- 


exp 


kA—0 kA,kB' 


II. CLICK MOMENTS FROM MEASURED CLICK STATISTICS 

Let us define the operators 


-.fri^fn^-. = ■,e-^AT(nA/N)^-lBT{nB/N).^ 

for 0 < li < N (i = A, B). Note that we use here the moments rhi instead oi fri = 1 — fhi, which were used in our 
earlier work [1]. The mathematical treatment in using fhi or -ki is identical. However, it is more convenient applying 
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nii for the present considerations. The matrix of joint click counting moments is 
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where the framed part includes the second order correlations It is denoted as the cross-correlation matrix, 

since detM^^’^) = — (:[AmA][Am_B]:)^ [!]• The single-mode reduced matrices of moments read 

as 
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Again, the framed inset includes solely the moments up to the second order, i.e.: and Those are defined 

as the variance matrices, because detM^^’°^ = (:[Atoa]^:) and detM*^°’^) = (:[AmB]^:). 

The moments may be retrieved via [1] 
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yielding an experimental mean value of 
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The standard estimate for the sampling error is 
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For 0 < Ka, Kb < N, the matrix of sampled mean values, ^ and the matrix of the errors, are 

given by the elements {'.rh^-l■) and a respectively, arranged according to Eqs. (5) and (6). 


III. DETERMINING CORRELATIONS 

For fully classical states, the matrices of moments are non-negative. There are genuine quantum correlations 
correlations of ATth order between the modes, if 

> 0, > 0, and ^ 0, (9) 

for 0 < K < N. This would mean that there are no non-classical correlation in the single mode marginals, but there 
are some nonclassical two-mode correlations. The non-negativity of a matrix M, or its violation, may be probed by 
its normalized eigenvectors v {Mv = ev, with v^v = 1 and e denoting the corresponding eigenvalue), since 

Mv = e. (10) 

If all eigenvalues of any matrix M are non-negative, Vn : Cn > 0, it holds that M > 0. Conversely, if at least one 
eigenvalue is negative, 3n : e„ < 0, we have a violation of the non-negativity M ^ 0. 

Experimentally, this may be obtained via the following steps: 















3 


1. determine the matrix of moments, i.e.: M ± AM; 

2. compute the normalized eigenvectors Vn of M to get 

3. perform a linear error propagation, which yields 

Ae„ = |n„|'l'AM|n„|, 

where \f \ = (|/i|, I/ 2 I,.. for any vector /; 

4. and, hnally, one gets e„ = e„ ± Ae„, or one may compute the signed signihcance = e„/Ae„. 

The minimal bound for the eigenvalues of classical states is Cd. = 0, cf. [1]. This means that the absolute value of the 
significance corresponds to the distance of e to this classical bound in units of error bars, |I]| = |e — Cd.|/Ae. The 
sign of S shows if we are consistent with this bound, S > 0, or clearly violate it, S < —3, meaning that the state is 
significantly nonclassical. Any significance level closer to zero cannot be well separated from the classical bound. 

Please also note that the presented approach has a favorable error propagation compared to the equivalent minor 
method in Ref. [1]. The conclusions from verihed negativities are identical with both techniques. Moreover, one should 
stress that the A^th order moments fully characterized the click counting statistics, which is of particular interest in 
case that second order moments yield inconclusive or insignificant results. 


IV. ALL RESULTS OF THE DATA ANALYSIS 


In this section, we present all information obtained during our data processing. We used the method as described 
above. In Fig. 1, the measured click statistics are given. In Tables I and II the eigenvalues of the single mode variances 
and the second order cross-correlations, respectively, are given. Similarly, Tables III and IV include the significance 
of the eigenvalues of the full matrix of moment for the bipartite and single mode case, respectively. 


50 //W 100 iMI 200 iMJ 



FIG. 1: The measured click counting statistics Ckj^,kB (for kA,kB = 0,..., 8) are shown at different pump powers 
{50 fiW, 100 fJ,W, 200 fiW, 300 /rW, 403 ^W}. A blank area depicts that no coincidence has been recorded. The total numbers of 
events are {41 750 447,41 755 519, 37 543 468,41 728 175,41 711 631}, being sorted according to the previous pump powers. The 
co,o coincidences are cut. They have the values {0.96,0.91,0.81,0.70,0.55}, respectively, yielding the major contribution to the 
whole statistics, i.e., more than 50%. The measurement time for each pump power was 10 minutes. 
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TABLE I: The eigenvalues (e^ , i = 1,2 and j = A, B) are shown for the single mode variance matrices and The 

minimal (boldface) eigenvalues are plotted in the main body. All eigenvalues are non-negative. 


pump power 

et 

ef 

e# 

ef 

50 pW 

1.9936867±6.5xl0"® 

1.9941778±6.3xl0"® 

(33.3±6.5)xl0"® (36.3±6.3)xl0"® 

100 mW 

1.985126±1.0xl0"® 

1.9863857±9.7xl0“® 

(9.2±1.0)xl0“® 

(97.3±9.8)xl0"® 

200 mW 

1.966550±1.6xl0"® 

1.968691±1.6xl0“® 

(27.0±1.6)xl0"® (28.2±1.6)xl0"® 

300 mW 

1.940763±2.0xl0"® 

1.945466±2.0xl0“® 

(64.0±2.1)xl0"® (63.0±2.0)xl0"® 

403 mW 

1.894942±2.7xl0"® 

1.897773±2.7xl0“® 

(165.7±2.9) X10"® (170.9±2.9) x 10“® 


TABLE II: The eigenvalues (ef^, i = 1,2,3) of the click correlations matrix are shown. The negative (boldface) 

eigenvalues are plotted. 


pump power 

^AB ^AB ^AB 

ei e 2 63 

50 pW 

100 mW 
200 mW 
300 mW 
403 mW 

3.976098±2.1xl0“® (148.4±9.8)xl0"® (-1.5±l.l)xl0"® 
3.944089±3.2xl0“® (39.3±1.5)xl0“® (-3.6±1.8)xl0"® 
3.873915±5.2xl0“® (111.0±2.5)xl0"® (-9.0±2.8)xl0"® 
3.780747±6.5xl0“® (243.5±3.2)xl0"® (-15.2±3.6)xl0“® 
3.608217±8.8xl0“® (604.6±4.7)xl0"® (-32.2±4.9)xl0“® 


TABLE III: The signed significances of eigenvalues (ef^, i = 1,, 10) of the full higher-order click correlations of the 
10 X 10 bipartite matrix of moments are given depending on the pump power. The third eigenvalues are negative with 

the highest signihcance (E 3 < —3) and, thus, clearly demonstrate the nonclassical character of the generated radiation field. 


pump 

power 

^AB \^AB \^AB ^AB ^AB v^AS v^AS ^AB ^AB 

2^1 Zj2 ^3 ^4 ^5 ^6 ^7 ^8 ^9 ^10 

50 pW 
100 mW 
200 mW 
300 mW 
403 mW 

e.OxlO"' 2.3x10^ -2.7x10° 2.2x10"^ -2.9x10"® -8.7xl0"'‘ 2.0x10“® 5.7x10“® 3.7x10“® -2.8x10“® 

3.9xl0'‘ 3.9x10^ -4.0x10° 5.0x10“® -6.0x10“® -1.3x10“® 7.2x10“® 2.9x10“® 9.2x10“® -6.9x10“® 

2.4xl0'‘ 6.5x10^ -5.9x10° 1.6x10“^ -1.5x10“® -3.0x10“® 2.1x10“^ -7.4x10“® 2.3x10“® 1.4x10“® 

l.OxlO"' 1.0x10® -7.5x10° 5.0x10“^ -3.2x10“® 2.0x10“® 1.9x10“® -2.5xl0“'‘ -8.5x10“® 1.5x10“® 

1.4xl0'‘ 1.7x10® -1.1x10® 1.7x10° -9.3x10“® 1.5x10“® -1.6x10“® -9.4x10“® 4.6x10“® 1.1x10“® 


TABLE IV: The signed significance of eigenvalues {el, i = 1,..., 5 and j = A, B) of the full higher-order moments are given 
for the single-mode reductions A and B. The single-mode matrices of moments and are of the dimensionality 

5x5. All significant values, E > 3, are non-negative, and most of the insignificant values are positive too. Note that any 
insignificant negative value, E < 0, is not inconsistent with the classical bound 0 and, thus, it cannot imply nonclassicality. 


pump power, mode 

V-? V-? V-? V-? 

Zji Zj2 2^3 2^4 2^5 

50 fiW, j = A 

8.7x10® 1.6x10® 2.3x10“® -2.4x10“® -7.7x10“®® 

50 mW, j = B 

9.1x10® 1.8x10® 7.5x10“® -4.9x10“® -8.4x10“®° 

100 jAN, j = A 

5.7x10® 2.8x10® 1.4x10“® -1.6x10“® 3.0x10“®° 

100 /iW, j — B 

5.9x10® 3.1x10® 2.7x10“® 2.4x10“® -6.3x10“® 

200 ^^W, j = A 

3.6x10® 5.0x10® 1.6x10“® -2.1x10“® -3.6x10“® 

200 fiW, j = B 

3.7x10® 5.3x10® 2.0x10“® 3.9x10“® -3.5x10“®" 

300 fiW, j = A 

2.8x10® 8.6x10® 8.8x10“® 5.5x10“® -4.1x10“° 

300 fiW, j ^ B 

2.9x10® 8.9x10® 7.7x10“® 1.4x10“® -9.2x10“®" 

403 /rW, j = A 

2.1x10® 1.5x10® 4.0x10“® 1.2x10“® -1.9x10“®" 

403 /iW, j — B 

2.1x10® 1.6x10® 4.0x10“® 8.8x10“® 1.4x10“° 
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V. FAKE EFFECTS 


The typically considered approximation is that the photoelectric detection model is valid for a large number of bins 
N and for low intensities. In other words, one assumes 


]kB 


Pk,,kB = \ -r(n.) [r(nB)] " ^-r(n^). 


kAl ks'- 

The moments of the photoelectric counting statistics PkA,kB &re given by 

kAl 


(:[r(n^r-[r(nBr-)= E E 


■) ~ ('kA,k 


knl 


kA>mA kB>mB 


{kA - mA)'- (ks - ms)! 


PkA,kB ■ 


( 11 ) 


( 12 ) 


In the simplest case of r(n) = n, these moments correspond to the normally ordered photon number moments, 
In order to prove that fake effects occur, we follow the same treatment as done for the click-counting 
moments. That is, we formulate the single mode 5x5 matrices of photoelectric moments of the marginals, i.e., 

and ((:[r(nB)]'""+“-(13) 

and compute the signed significances and S® of the minimal eigenvalue. If the approximation (II) was true, then 
the marginal thermal states must not exhibit significant negativities. The corresponding falsification with our data is 
given in Table V. This simply means that the approximation in Eq. (II) is not justified for the used detection scheme. 


TABLE V: The minimal signed significances {j = A,B) of the eighth-order photoelectric moments are given for the 
single-mode reductions A and B. The negativities demonstrate that the photoelectric detection model is not adequate for our 
detection scheme. They tend to increase with increasing pump powers. Apparently the fake effects are highly sensitive to the 
modes - label A corresponds to the earlier time bins and B to the later ones - which might be due to an overestimation of 
afterpulses using the wrong measurement statistics. 


pump power 

E'^ E® 

50 mW 

-2.6 -2.4 

100 pW 

-3.1 -1.4 

200 pW 

-2.6 -2.2 

300 AtW 

-6.6 -3.1 

403 pW 

-12 -13 


VI. THE THEORETICAL MODEL AND THE INFLUENCE OF IMPERFECTIONS 


In order to compare the measured results with a theoretical prediction, we outline the corresponding theoretical 
model. The target state is a two-mode squeezed-vacuum state (^ > 0), 


i(-\ (tanh^)’’ 1 ^ 1 ^ 

= E _uc mA\n)B, 


n—0 


cosh^ 


(14) 


where the squeezing parameter is a function of the pump power: ^ = ^o^/P. The relation to the energy per pulse is 
T-pump = P/70 kHz, for the given repetition rate. The detector response T in a linear form is 


r{n/N) = rjfi/N + v. 

The resulting parameters are fitted: 

^0 = 0.087 (/xW)-i/2^ p = 0.096, and = 0.51. 

The relation between the squeezing parameter and the mean photon number is 

n = (^|nA + nsIC) = 2 sinh^ ^ or ^ = arsinh • 


(15) 

(16) 

(17) 
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In Table VI, the pump power and the resulting energies i?pump) squeezing parameters mean photon numbers n are 
listed. The squeezed quadrature of this state is X = xa — xb = [a + a’^] — [6 + &’*’] yielding to (AV)^ = (^|[AA]^|^) = 
2e~^^ and a noise suppression of 

® 

However, the measurement setting under study cannot be used to infer phase information, because it is designed to 
perform photon correlations measurements. Hence, this theoretical relation does not show the true squeezing level, 
for example, phase randomizations can signihcantly reduce squeezing levels and are not accessible with our scheme. 
A similar conclusion holds true for entanglement. 

TABLE VI: For the pump powers P, the following quantities are listed: the energy, Epump = P/70 kHz, the squeezing 
parameter, f = 0.087 and the mean photon number, h = 2(sinh[0.087 (/rW)“^^^\/P])^. 


p 

50 mW 100 mW 200 300/rW 403/rW 

-^/pump 

0.7 nJ 

1.4 nJ 

2.9 nJ 

4.3 nJ 

5.8 nJ 


0.6 

0.9 

1.2 

1.5 

1.7 

n 

0.9 

1.9 

4.9 

9.2 

15.5 


For a comparison with the measured data, the moments {irh^rh^^:) 


can be calculated numerically or analytically: 




^-v(kA + kB) 

cosh^ ^ — [1 — -qkA/N] [\ — rjkB/N] sinh^ ^ 


(19) 


It can be directly seen that the dark count rate is just a scaling of the moments. For a given v, it scales the eigenvalues 
but does not influence the significance of determined negativities, because the errors scale in the same way. 

In Fig. 2, we plot the minimal eigenvalues of the theoretically obtained second order matrix of click moments. Due 
to symmetry, we have equal eigenvalues for the single-mode subsystems A and B. For very small and very large 
mean photon numbers h, the minimal eigenvalues converge to 0. For small n this is easy to see because we approach 
the vacuum state. For large n, the mean photon number increases so that saturation effects of the APDs diminish 
correlations, which is correctly displayed by the theory via the decay of correlations. 



FIG. 2: The minimal eigenvalues ca = es of (left) and cab of (right) are plotted depending on mean 

photon number n and the quantum efficiency rj (fixed dark count rate v = 0.51; for a better visualization, we plot —cab). The 
single-mode eigenvalues ca = cb are positive for any value of h > 0 and rj > 0, while the bipartite case shows negativities for 
all these parameters, —cab > 0. 


[1] J. Sperling, W. Vogel, and G. S. Agarwal, Phys. Rev. A 88, 043821 (2013). 













13 



